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Abstract 

We study the gravitational corrections to the F-term in four-dimensional M = 1 
U{N) gauge theories with flavors, using the Dijkgraaf-Vafa theory. We derive a 
compact formula for the annulus contribution in terms of the prime form on the 
matrix model curve. Remarkably, the full Tri? A R correction can be reproduced 
as a special momentum sector of a single c = 1 CFT correlator, which closely 
resembles that in the bosonization of fermions on Riemann surfaces. The N = 2 
limit of the torus contribution agrees with the multi-instanton calculations as well 
as the topological A-model result. The planar contributions, on the other hand, 
have no counterpart in the topological gauge theories, and we speculate about the 
origin of these terms. 



^hfujiOpost .kek.jp 
^mizoguchOpost .kek.jp 



1 Introduction 



Since the discovery by Dijkgraaf and Vafa, we have recognized that non-perturbative 
aspects of four-dimensional supersymmetric gauge theories can be studied via matrix 
models. In this framework the effective superpotential for M = 1 supersymmetric gauge 
theories can be determined as the large free energy of a matrix model and by 
minimizing it the non-perturbative vacua and their phase structures can be investigated 
P] . This part of the proposal has been elegantly proven by using the chiral ring relations 
of = 1 supersymmetry in the generalized Konishi anomaly equations OH]. 

On the other hand, the non-planar diagrams have been shown to correspond to the 
gravitational corrections to the gauge theory El 13 El El CDI, in particular, the genus 
one free energy of the matrix model computes the gravitational F-term 



Recently, there has been much progress in the analysis of A/" = 2 gravitational F-terms in 
terms of the multi-instanton calculations ^T] and geometric engineering IT^ ITlj . To 
compare them with the matrix model results, it is essential to know the precise form of the 
gravitational F-terms computed using the matrix model. For the pure gauge theory case, 
it has been checked that the torus free energy J^i coincides with the topological partition 
function J3] of the Donaldson- Witten theory, and hence the M = 2 gravitational F-terms 
iniin]. The planar free energy is also known to contribute to the gravitational corrections 



This framework can be extended to A/" = 1 supersymmetric gauge theories with Nf 
flavors in the fundamental representation [T71 ITHl IT^ . The dual matrix model consists of 
a bosonic matrix and vectors which correspond to the adjoint scalar cf) and Nf flavor fields 
, qi (/ = 1, ■ ■ ■ , Nf) of the gauge theory, respectively. In the vector coupled matrix 
model, there are two kinds loops with N and Nj indices; the former are summed up but 
the latter remain as boundaries of Feynman diagrams. It looks as if the model has not 
only the closed string but also the open string sectors [201 • In order to evaluate E? terms 
for this gauge theory, it is necessary to consider both the torus and annulus contributions. 

In this paper, we evaluate the correction terms for A/" = 1 U{N) gauge theories 
with Nf flavors. Using the ordinary large analysis [21], we derive a compact formula 
for the annulus contribution to the correction in terms of the prime form on the matrix 




(1.1) 
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model curve. Remarkably, the full TiR A R correction containing the torus as well as all 
the planar contributions is reproduced as a special momentum sector of a single c = 1 
CFT correlator, which closely resembles that in the bosonization of fermions on Riemann 
surfaces. It is in accord with the recent observation in topological string theory that the 
non-compact B-branes correspond to fermions in a chiral boson theory on a hyper-elliptic 
curve 1221 123- 

The plan of this paper is as follows: In section 2 we briefly review the Dijkgraaf-Vafa 
theory including a string theory derivation of gravitational corrections. In section 3 we 
compute the planar gravitational corrections for the A/" = 1 U{N) gauge theories. In 
particular, we derive in section 3.4 a compact formula for the annulus contributions in 
terms of the prime form on the matrix model curve. Section 4 is devoted to some examples, 
in which our formula is confirmed explicitly. In section 5 we show that the full TiR A R 
gravitational correction, including both the torus and planar contributions, is reproduced 
as a single chiral correlator of a c = 1 conformal field theory, which closely resembles that 
in the bosonization of fermions on Riemann surfaces. In section 6 we consider the Af = 2 
limit, and speculate about the origin of the planar contributions. Finally, we conclude our 
paper with a summary and outlook on future work in section 7. The appendices contain 
some technical details which we need in the text. 

2 Gravitational F-terms from Matrix Models 

2.1 The Dijkgraaf-Vafa Theory with Flavors 

The original proposal of Dijkgraaf and Vafa[T] was summarized as follows : 

1. The low energy effective superpotential of A/" = 1 gauge theories can be computed 
by summing over the planar diagrams of the matrix model with the same tree-level 
potential. 

2. The non-planar diagrams compute the gravitational F-terms for these theories. 

The first statement was proven in [3], and was generalized to the cases with flavors in 
[THj . The latter part was also supported by many arguments [Hlini, and was explicitly 
confirmed in the pure gauge theory cases [31 IH] . 
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Let us consider the = 1 U{N) gauge theory coupled to matter superfields qj, qi 
(/ = 1, ■ ■ ■ , Nf) in the ( ant i-) fundamental representation with the superpotential J7j 

V{(j),qi) = tTWi(l)) + J2qi{(f)-mi)qi, 

1=1 

W{^) = E^r- (2.1) 
p=l p 

4> is the adjoint chiral superfield. In the classical vacuum the gauge group is broken as 

n n 

U{N) ^ UUiNi), Y.Ni = N. (2.2) 

1=1 1=1 

The claim of [T is that the non-perturbative vacuum structure of this theory can be 
analyzed by a vector coupled matrix model with action 

5matnx(<f , Qi, Q i) = Tt^^^W{^) + ^ Qj{^ - mi)Qj (2.3) 

1=1 

where $ is an iV x iV hermitian matrix and Qj, Qi are iV-component vectors. The vector- 
matrix coupling Qi^Qi leads to Feynman diagrams with boundaries in various topologies. 




Figure 1: Feynman diagrams with boundaries. 

In the low energy effective theory, the glueball superfields Si = -^^tr su{Ni)WaW° 
[i = l,...,n) play the role of the fundamental fields. According to the proposal, Si 
is identified with the 't Hooft coupling QsNi for each i, where Qs is the matrix model 
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coupling constant and Ni is the number of eigenvalues distributed on the i-th cut. Under 
this identification, the effective superpotential for the gauge theory is given by 

W^siS.) = ± (n^^^Q^ + rs}j + 4'\S.) (2.4) 

where ^o°^(S'j) and J-'Q^\Si) are given by the large N sphere and disk free energies of the 
matrix model, respectively, r is the bare coupling.^ By extremizing this superpotential, 
one can analyze the non-perturbative vacuum structure of the gauge theory j2]. 

2.2 String Theory Derivation of Gravitational F-terms 

To extend the analysis to the study of the gravitational F-terms, it is useful to consider 
how they arise in string theory. The gauge theory setup above can be realized in type IIB 
string theory on a Calabi-Yau three-fold [211 123 

W'izf + y'^ + + = 0, {z,y,v,w) eC^ (2.5) 

with n singular points. By blowing them up, we obtain a smooth geometry with excep- 
tional 2-cycles Ci {i = l,---,n). We then consider D5-branes wrapped around Cj's 
and Nf D5-branes around non-compact 2-cycles C/ (/ = 1, ■ ■ ■ , Nf). The Af = 1 gauge 
theory above is realized on the space-time filling world-volume of the D5-branes. 

To evaluate the F-term corrections it is advantageous to utilize the hybrid formalism 
j26t I7j. The string Lagrangian density for the four- dimensional part is given by 

c = ^dx^'Bx^ + pje'' + p^Be^ + pje" + paOe^, (2.6) 

where p, p, 9 and 9 are the fermionic fields. Inserting two gluino vertex operators §^ W^Pa 
on the boundaries 7 of the world-sheet, the stringy corrections lead to the F-terms con- 
taining the glueball superfields. 

The string world-sheet can end either on compact or Nf non-compact branes in 
the Calabi-Yau direction. Resorting to the fermion zero-mode arguments |771 1^ . we 
can conclude that only the planar string world-sheets with at most one boundary ending 
on the non-compact branes can contribute to the effective superpotential jT]. A simple 

^The definition of Wcfi{Si) H2.4|l has an ambiguity in tlic linear terms in Si depending on tlie clioices 
of the integration constants in the free energy (See e.g. jl8|.). 
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Figure 2: The brane configuration in the resolved geometry. 

combinatorial argument then yields the effective superpotential ()2.4j) where the sphere 
and disk free energies of string theory are identified with those of the matrix model. 

The gravitational F-term corrections are given by inserting RR vertex operators on the 
bulk of the string world-sheet. A candidate for such an operator is the self-dual gravitino 
vertex operator [3 |H] 

J G'^^^I^Po.Xp^dX.^+p^X^^dX.^y^^ + J G'^f-'p^ppid, - e,). (2.7) 

Due to the chiral ring relations 

{W^,Wp} = 2Gc,p^W\ G^ = Q (mod 5), (2.8) 

at most two Gap-^ insertions can contribute to the gravitational F-term corrections. Such 
G^ terms are given by the torus (Figure El) as well as the planar world-sheet corrections 
(Figure 1^. In particular, there are three types of planar world-sheets depending on which 
(compact or non-compact) branes the string ends on. 
The gravitational F-term containing G^ takes the form 

j d*xd^ej^^=o{Si)G^, (2.9) 

where J-'^=o{Si) (= the planar + torus gravitational corrections) is the contribution from 
the X = string world- sheets.^ If we further consider a self-dual graviphoton background 
^Here only the empty loops are understood as boundaries in the definition of the Euler number x- 
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Figure 3: Gravitational corrections from the torus. 



Fapi then we also have gravitational F-terms from the graviphoton fields. Such F-terms 
are obtained by decomposing M = 2 gravitational F-terms into the M = 1 F-terms. It 
was shown in [2j that the superstring computation of these Feynman diagrams is identical 
to that in the field theory limit, and the combinatorial factors can be calculated in the 
associated matrix model. 

In the following sections, we will derive a formula for the free energy T-^=q{S,i). 

3 Planar Gravitational Corrections 

3.1 Gravitational Corrections from Sphere and Disk Diagrams 

There are three kinds of planar diagrams which contribute to the Tri? A R correction 
(Figure E}. The first and second types of contributions are determined by the sphere and 
the disk free energies through the following combinatorial arguments: For the unbroken 
gauge group case (one-cut case), the coefficient of planar gravitational corrections yields 

oo A/-2 f)2 -77(0) ( c\ 00 

h=l ^ h=l 
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Sphere with 2 index loops Disk with 1 index loop Annulus 



Figure 4: Three types of planar string world- sheets, 
for the first type, and 

E^c^:F^'\^^suiN)irs'^-' = N^^^^ :Ft,'\s) = j:^l^is'^ (3.2) 

h=l h=l 

for the second type. J-'Ifl and J-'^^l are the symmetric factors of Feynman diagrams with 
h boundaries for the sphere and disk topologies, respectively. Generalizing to the cases of 
arbitrary breaking patterns, we have 

1 " d'^T^^^ 

- ^ ^i-^j an an (^'-'^ ^ sphere with 2 index loops), (3.3) 



2,^, ^dS,dS, 
' dS, 



Nj ° (for a disk with 1 index loop). (3.4) 



i=l 

The third type of gravitational corrections come from annulus diagrams, which we will 
compute below using the familiar large N technique in matrix models j2l] . 

3.2 Planar Diagrams with Boundaries 

Let $ be an iV X iV hermitian matrix, and {I = 1^ . . . ^ Nf) be iV-dimensional complex 
vectors. We consider the partition function 



vol U{N) 
X exp 



1 

9s \ 



1=1 

Nf 



1=1 



(3.5) 



with a polynomial tree level potential W{^). It allows a topological expansion j2H] 

T = >^>^gf+'=^f, (3.6) 



g=(}k=0 



where JF^'^) is the contributions from the connected genus g diagrams with k boundaries. 
J-'q^^ also includes the non-perturbative piece coming from the gauge volume pUj . 

(2) ~ 

To evaluate the annulus amplitude JFq , we carry out the dQjdQi integrations first. 
This yields 



oc 



(i$ exp 



Nf 



Tr 

gs 



1=1 



(3.7) 



up to an irrelevant multiplicative constant which is independent of $, iy($) or m/'s. 
This integrations organize a resummation of diagrams. The gs factor of the log potential 
indicates that there is a matter (= QQ ) loop at each log vertex. Thus, instead of Z, we 
consider the following double expansion 



Z{gs 



^{gs 



1 



Nf 



Tr 

gs 



(i$exp 
exp(^7^J^(^„e)), 

oo 

9=0 

oo 

k=0 



W {<!>) + eJ2\og{<!> - mj) 



(3.8) 

(3.9) 
(3.10) 



and set e = gg- Then J^f'^ and JFq^-* are of same g^ order in the ordinary g^ expansion of 
the free energy; the former is the torus, whereas the latter is the annulus amplitude. 

One of the aims of this paper is to determine the precise form of JFq for general 
n and for a general polynomial potential W{z), in terms of the language of Riemann 
surfaces that the matrix model defines. JFq^^ [k = 1,2, . . .) can be extracted from Z{gs, e) 
as follows: Writing the expectation value of any function /($) of $ with respect to the 
measure (j3.8j) as (/($))^^ ^, we have 

(fc) ^ 1 g"^o(e) 
k\ dt^ 



e=0 
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1 ,„„s'^te.^) 



e=0 



— lim 



k~l 



k\ de 



k-l 



9s(TTEfli log{^- mi) 



(3.11) 



e=0 



Keeping S := QsN fixed, — > implies N —>■ oo. In this limit the expectation value can 
be written as 



lim i(Tr/(<l>))^ 



rfAp(A,e)/(A) 



for any /($), where p(A, e) is the eigenvalue density function of $ normalized as 



dXp{X, e) 



Therefore, expanding p(A, e) as 



p(A,e) = E^V'HA), 

fc=0 



(3.12) 



(3.13) 



(3.14) 



we find 







(fc) 



1 (9' 



— -S / rfAp(A, e) log(A - m/ 

: \ J-oo 



/ dAp('=-i)(A)^log(A-m,), 

K J— oo r 1 



e=0 



1=1 



in particular, the disk [TTl HH] 



oo 

^(S'^ = rfAp(°)(A)EMA-m,). 

J-00 



and the annulus 



4'^ = -f rrfAp«(A)EMA-m,). 

/ J~-oo 



(3.15) 



(3.16) 



(3.17) 



p('S"i)(^;\) can be computed by means of the standard large N technique for evaluating the 
planar diagrams fH] . 

To be rigorous, one would need to show that the Qs ~^ limit and the operation of ^ 
commute. In Appendix A we will give a different derivation of the formula ()3.17|) for an 
arbitrary polynomial potential Vr($) using Riemann's bilinear identity, thereby providing 
an alternative proof of it. 
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3.3 The Large N Technique 

We will now briefly review the large (= N) technique developed in (21] to fix our 
notation. (See also As we introduced above, p{X,gs) is the density function of the 

eigenvalue A of $ normalized as ()3.13p . The nonzero support of p{X,gs) (as a function of 
A) consists of several disjoint intervals on the real A axis. These intervals are the branch 
cuts of the resolvent uj{z, gs) defined by 

uj{z,gs) := / dX — (3.18) 

J — oo Z — A 

on the complex z plane. ()3.18|1 is equivalent to 

ij{X + iO,gs) -uj{X-iO,gs) = -2Trip{X,gs) (3.19) 

in the boundary value representation. If the number of cuts is n, uj{z,gs)dz is a mero- 
morphic differential on a hyper-elliptic curve of genus n — 1 with a parameter gg. 
Let Sn-i be a family of genus n — 1 hyper-elliptic curves defined by the equation 

2n 

y^ = l[{z-a,{gs)) {i = l,...,2n) (3.20) 

i=l 



with a parameter gs- We denote by yH^Zii^ — ai{gs)) a branch of y such that y > when 
2; e R, z ^ 00. We call it the first sheet, and the other the second sheet. 
The n cut solution to the large saddle point equation 



Lj(z,gs) = - — — > f aw , (3.211 



defines a meromorphic differential on Sn-i- The contour surrounds all the cuts but not 
the points w = z, mi {1 = 1,..., Nf). The definitions of the contours are summarized in 
Figure |31 The positions of the 2n branch points ai{gs) {i = 1, . . . , 2n) are determined by 
the 2n conditions 



N 




X\W'{w) + Y.lil^r 



J2f \dw \ ' ' ^-^u,-m,^ =AmS5k,n {k = Q,...,n) (3.22) 



and 



S- 

- dw u{w,gs) =2Tci-^ {j = l,...,n). (3.23) 

A, O 
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Figure 5: The contours. Those on the second sheet are shown in broken hues. is 
defined to go along \z\ = |Ao| counter-clockwise on the 2;-plane. 
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The first n + 1 conditions (j3.22j) arise from tlie requirement tliat uj{z,gs) beliaves like ^ 
as 2 ^ oo on tlie first sheet, whereas the latter are the conditions for a given set of filling 
fractions. Only n — 1 of ()3.23|) are independent. 
The large N free energy for the potential 

Nf 

W{^,gs) := W{^)+gsJ2^og{^-mj) (3.24) 

1=1 

is given by 

/oo 
dXp{X,g,)W{X,gs) 
-oo 

/oo roo W \'| 

dX / dX'piX, g,)p{X', gs) log (3.25) 
-oo J —oo 1\ 

A is an arbitrary integration constant, which may be regarded as the physical scale pa- 
rameter of the corresponding gauge theory up to a potential dependent constant. 

3.4 Annulus Diagrams 

We would now like to write ()3.17|) in terms of the language of Riemann surfaces. ^^(A, gs) 
is expanded as 

oo 
k=0 

The differential uj^^\z)dz must have zero A-periods because the right hand side of (|3.23p 
does not depend on gg. 

Let us analyze the singularities of the differential u^^\z)dz. Since the discontinuity of 
uj{z,gs) is the eigenvalue density, it satisfies 




(3.27) 



u!{z,gs) also satisfies 

dzu;{z,g,) = l (3.28) 

because u{z,gs) goes to ^ at infinity on the first sheet. ()3.27p()3.28|) imply that u{z,gs)dz 
has a pole of first order at 2; = oo on the first sheet with residue —1, but otherwise it is 
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regular everywhere else on the first sheet. Therefore we conclude that u}^^\z)dz is regular 
everywhere on the first sheet. On the other hand, uj{z,gs) can be written in the form 

^iz,9s) = ^(w'{z) + j:-^-y{z,9. 



y{z,gs) = — If -\f \dw 




V^^- /=i^^V {z-w)^Utii^-a^igs)) 

(3.29) 

We similarly expand 

oo 

y{z,gs) = j:9sy^'\^)- (3-30) 

A:=0 

Note that the difference of the contours of y{z, gg) and uj{z, g^)] thanks to this change, the 
contour integral of y{z,gs) gives some rational function of z.^ This implies that y{z,gs) 
changes its sign under the interchange of the first and second sheets. 

Since the second term of u!{z,gs) has a first order singularity at each z = mi (/ = 
1, . . . , Nf) with residue |^ and at 2; = oo with residue —^^^ on both sheets, the y{z, 
term must cancel them on the first sheet, and hence in turn doubles the residues on 
the second sheet. To summarize, denoting coordinates on the second sheet with tildes, 
uj^^\z)dz has first order singularities at 

• z = rhj with residue ^ for all / = 1, . . . , Nf, and 

• z = 60 with residue — 

otherwise holomorphic everywhere else on Sn-i- Thus we have shown that uj^^\z)dz is 
an abelian differential of the third kind with zero A-periods and is given by 



UJ 



1 



Z)dz = -J2^nii-6b, (3.31) 

1=1 



where, following the notation used in ^22], uja-b{z) denotes a zero A-period abelian differ- 
ential which has simple poles at z = a,b with residue ±1, respectively. 

One of the nice things for this kind of differentials is that their integrals are compactly 
described in terms of the prime form E{z,w) on S^-i- The prime form E{z,w) is known 

■^If gs — 0, it is a polynomial of z. If, in particular, n saturates the degree of W'{z) (the maximally 
broken case), it is a constant up to which y^^'>{z) coincides with y{z) (|3.20() . 
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to be the unique bi-holomorphic half differential on a compact Riemann surface such that 
E{z, w) = if and only if z = w. One of the basic properties of the prime form is 

E{z, a) 



dz log ■ 



(3.32) 



E{z,h) 

Using this formula, we finally obtain the third type of contributions to the gravitational 
correction^ 

_ ^ / 





-77(2) 

■J n 



4m ^ 

1 ^ , E(mj,6o)E(oo,mj) 

= — y log z — 

2 j'j^^ E{mi, mj)E{oo, oo) 
The full planar gravitational correction JF(p^) is thus 



(3.33) 



(pi) 



with 



r 



in _q2 -r-(O) 

(1) 



1 n q2 -r-(O) n /^t(1) 



lim 

|A„Hoo 



E N.N, 



^,(z)+log^ 



(3.34) 



(3.35) 



EN. 

1=1 



n " f 

J2nJ uj,{z)J2\ogiz- 
,fc=l "^^fe 1=1 

i=i 1=1 \ -''^i 



mi) 



lim 

lAoHoo 



loef-Ar 



(3.36) 



and Tq'^ given in ()3.33|) . In deriving the above, we have used the special geometry relation 

PHHI 



dSi 



lim - 



|Ao|-*oo \2 Jb. 



dz y^°\z) -W{Ao) + 2S log 



-A 

A 



(3.37) 



^Of course, the computation of the annulus amplitude J^q^' itself can also be (and has been) done by 
other means (See e.g. [331 1^.). What is new here is that we have expressed it in terms of a special form 
(that is, the prime form) on the Riemann surface, and we cannot compare it with the CFT result until 
we do so. We thank H. Kawai for discussions on this point. 
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as well as the variation formula 



dz = —4:mu;i{z){i = l,...,n), 



where cui {i — 1, . . . ,n) are related with the holomorphic differentials cui {i — 1, 
on Sn-^^ through the relations 

1 

= - 7r^<^oo-6b (^ = 1), 
2m 

^ _ 1 
2m 

They are normalized so that 

Uj(z) = 5ij = 1), 
u}j{z) = 6ij = l,...,n). 



A, 



See Appendix B for details. 

4 Examples 

4.1 One-cut Solutions for Quadratic Potential with Nf = 

The first example is the 'conifold' case 

with a single matter field. The resolvent u!{z,gs) is 

9s 



- 2tz + 



z — m 



+^(z-c{gsW-^i{gsy{2t 



^J {m - c{gs)f - liigsY 



where c{gs) and fj.{gs) are related with the positions of the end points by 

ai(9s) = c{gs) + i^{gs), 
a2{gs) = c{gs)-n{gs) 
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and satisfy (jlO^ and 



= 2tcigs) + ^= (4.4) 
^(m-c(^,))2-/i(^,)2 

25 = -t^,ig,r + g.\l- , |. (4.5) 

V("^-c(^s))2-/i(^,)^ 



They are solved as a power series 



c(^7.) = ; ; . ^.+ i;; ',:"^ 9s + o{g!), (4.6) 



y/m'^ ~ fi"^ — m —{m^ — 2/i^ — rriy/m'^ — fi'^Y + 5fi^ a ^-i/ s 

^^^^^ = ^+^w^^^^ + 4tw{m^-^^^r 



(4.7) 

with /i := The integral for the genus zero free energy jFg ()3.25|) can be completely 

performed in this case and is given in Appendix C. JFq^-* (the disk amplitude) and JFq^^ 
(the annulus amplitude) are obtained as coefficients of the Taylor series of J-'o{gs)- After 
some calculations using ()4.6|1 and ()4.7j) we find 

= -g(log '"+^f^ + " ^ -1). (4.8) 

\ 2 m + v?^ /U"^ 2y 

= ilogI^±^^ = (4.9) 



2 2Vm2 - /i2 V t 

agrees with the calculations in |19j . 
Let us show that ()4.9|) can be written in terms of prime forms. Let w,w' G CU{cx3} = 
P-*" be coordinates of some two points in the natural coordinate system, then the prime 
form is simply given by 

w — w' 

E{w,w') = Z_ "_ ■ (4.10) 



dwVdw' 

is also realized as a two-sheeted Riemann surface with a single branch cut. We take 
two end points at z = ±/i (/i G R, /i > 0), where z is the coordinate of the two-sheeted 
Riemann surface. The relation between these two coordinate systems is 

/U^ Z + y/z'^ - /i2 / . , , X 
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This map is two to one for generic except at the two end points of the cut where it is 
one to one. 

The cut [ — /i,/i] on the two-sheeted z plane gets mapped to a circle |w| = | on the 
w plane. Defining the first and second sheets as before, the region outside this circle 
corresponds to the first sheet, while inside the second sheet. It is easy to see that 



E{zrz) 



Aw 



dw (z) ^ dw (z) 



E{z, oo) 



dw{z)^dw{z) 
w-0 
^dw{z)dw{(X)) 



E(oo, z) 



^Jdw{z)dw{oo) 

oo-f 
^dw{oo)dw{z) 



E{oc,6o) = ^ =, (4.12) 

Jdw{oo)dw{oo) 

where, as before, we have added tildes to the z coordinates for the points on the second 
sheet. Thus one may write (j4.9|) as 

= Log ^K^W--'") (4.13) 

2 E[m,'m)E [00,00) 

This agrees with our general formula (j3.33|) . 

4.2 Two-cut Solutions for a Quartic Potential with Nj = 2 

The next example is an even quartic superpotential with two flavors with masses 

mi = —1712 '■= TTi-, (4-14) 
in which we find a symmetric 2-cut solution. The potential is 

W{^) = - u^^. (4.15) 
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Without matter, a large symmetric solution was obtained long time ago by ^H] in 
terms of elementary functions. Despite the symmetric potential, asymmetric filling is also 
possible and the solution is written using elhptic functions in general [SSIEZI- 
After integrating out the matter fields, we have 

W{\ gs) = -tX' - uX^ + gs \og{X^ - m"). 

With a symmetric ansatz, we find the resolvent 



2S'u;(z, gs) = -2tz - Auz^ 



9s 



9s 



z — m z + m 
4u 



'29s 



[z'^ — m?)J {w? — a'^){m? — 6^) 



,(4.16) 



where two cuts are created at [—6, —a] and [a, 6]. The density function of eigenvalues 
yields 



p{\9s 



2nS 



XJ {X^ - o?){h^ - X') 



-Au + 



2gs 



■ (4.17) 



Being symmetric, the condition (j3.23|) is automatically satisfied, while the asymptotic 
conditions ()3.22|) yield 

9s 



= 2mc + t + 



where parameters c, /x and M are defined by 



(C - M)2 - /i2 

M-c 



,y(M-c)2-^2 



(4.18) 
(4.19) 



a + h b — a , 

c = , u = , M = m . 

2 ' ^ 2 ' 

From these conditions we find an iterative solution 

2g'iu{2t + AuM - J{t + 2uMy + AuS) 



t 



9s 



+ 



/i 



2^ ^{t + 2uMY + AuS ' [it + 2uM)^ + AuS]^ 

t + 2uM \ 



+ 0{g, 



s 



{t + 2uMf + AuS^ 
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9t ru 



+ 



t + 2uM 



it + 2uMf + AuS AuS 



{t + 2uMf + 4:uS ^ 



+ 



2(t + 2uM){2t + AuM - J{t + 2uMf + AuS) 



+ 0{gl). 



(4.20) 



The multi-cut large N free energy is given in Appendix A, and is in the present case 
^0 = -S J dXp{X,g.,) (^liy(A,^,)-|log '^'~/°' ^ -|py(Ao,^.), (4.21) 

where jZbdXp{\gs) = Jl^dXp{X,gs) = 1/2. Plugging (HTTj) KM into and picking 

up the 0{g'^) terms, we obtain the free energy of annulus diagrams as before. The result 
is 



f2^ 1 t + 2uM + J{t + 2uMy + AuS 
-^0 =2 



2^{t + 2mM)2 + AuS 
Let us compare this with ()3.33|) . The prime form is defined on the curve 

y' = {z' - a'){/ - b% 

2n V M 2u ^ u 

In the present case ()3.33p reads 



(4.22) 



(4.23) 



1, E(m,m)E(m, —m)E(—m,m)E(—m, —m) ^ f f°°\ , x 

(4.24) 

The relevant meromorphic 1-form with zero A-periods is 

_ 1 / y{m) 1 ?/(-m) 1 1 I \ zdz 

l^m-6b + ^(-rh)-ao — " ^ — T^; 1 ; ; "2: -— . 



2 \ y(z) z — m y{z) z + m z — m z + m 



(4.25) 



By changing the integration variable ^ to x = 2;^, the above integration reduces to a 
simple form^ 



2 Jm \x — M y{x) X — M y(x) 



^Identifying z ^ —z, we obtain the same curve and singular points as those in the conifold case. 
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1 t + 2uM + J{t + 2uMy + AuS 
- log , ^ . (4.26) 

2 \ 2J{t + 2uMf + AuS 



This result exactly coincides with the matrix model calculation. 

4.3 Two-cut Solutions for a Cubic Potential with Nf = 1 
4.3.1 Perturbative Computations in Gauged Matrix Models 

The final example is the perturbative computations of the two-cut free energy for a cubic 
potential jSH] 

= c/Tr + — + g($ - m)g. (4.27) 



.3 2 

In this two-cut case we need to consider the fiuctuation around the vacuum 

^ = { 1° U f ) ■ (4-28) 

where ai and 02 are classical roots of W'{z) = 0. Around this vacuum the original gauge 
symmetry of matrix model 

$^f/.$.f/-i^ UeU{N) (4.29) 

reduces to U{N) U{Ni) x U{N2), and the matrix model action is given by 

W,U^) = gTr (1$?, + 1$?^) + gTr (^Ul, - 1$^,) . (4.30) 

Since the off-diagonal block of the matrix does not appear in the action, a convenient 
gauge choice is 

^^J=0, t^J. (4.31) 

In this gauge the coupling to vectors Q,Q is 

W^mattcr = Qi(*ii " m)Qi + 0^(^22 -m- A)Q2. (4.32) 

The gauge fixing requires the introduction of the ghost matrices B and C with action 

W^ghost(5,C) = ATr(52iCi2) - ATr(5i2C2i) 
+Tr(52i<l>nCi2 + C2i^nBi2) 

+ Tr (5i2$22C21 + Ci2$2252l) . (4.33) 
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4) Q B 




4> Q 



Figure 6: Vertices in the gauged matrix model. 



Figure 7: The dump-bell diagram. 
From the matrix model action 

W{<!>,B,C) = W^trce($)+W^ghost(5,C), (4.34) 

we can read off the propagators 

(^ll^ll) = ($22*22) = 

(-B12C21) = — (-B21C12) = — , 

9^ 

{QiQi) = --, {Q2Q2) = (4.35) 

m m + A 

For $^ and B^C vertices we assign weight g, and for Q^Q vertex we assign weight 1 
(Figure El). For a ghost loop, we add an extra factor of —1. 

Having found the Feynman rules, we can calculate the annulus contributions per- 
turbatively. The annulus diagrams are drawn in Figures [71121 up to three loops,. The 
contributions to annulus free energy from the dump-bell diagram Figure [7| is 

5*1 S2 



2gAm? 2^A(m + A)2' 



(4.36) 
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From the three-loop diagrams, we obtain the terms involving the Sf, 5*1 and 5*1 5*2 . 
The diagrams in Figure |H1 lead to the Sf and S"! terms. We obtain the Sf terms 



1 



1 



+ 



1 



+ 



1 



+ 



1 



+ 



1 



g^A^m^ g'^A^m^ Ag'^A'^m* g'^A'^m* g'^A'^m? g'^A'^m? 



Si (4.37) 



The 5*1 terms can be obtained by replacements ^ ^2, A ^ — A and m — * m + A 



+ 



+ 



+ 



^2A3(m + A)3 ^W(m + A)3 Ag'^A'^{m + AY ^W(m + A)4 
1 1 



+ 



+ 



Si 



^2A4(m + A)2 g^A\m + Ay\ 
Finally, the diagrams in Figure El lead to the S1S2 terms. 



g'^A^\m? {m + Af) g^A^ \m? {m + Af 
2/1 1 \ 2 



^f^A-^ y rrfi {m + A)^ J g'^ A^m{m + A) 
4.3.2 Comparison with the Annulus Formula 



S'i5'2. 



(4.38) 



(4.39) 



To compare these diagrammatic computation with the annulus formula ()3.33j) . let us 
expand the latter in powers of 5*1 and 6*2. For the two-cut solution the prime form is 
defined on the curve 



W{z) 



-W'iz) + h{z) = {z- z^){z - ^2)(^ - z^){z - ^4), 



/I 3 ^2 

nr + 2^ 



(4.40) 



There are two 1-cycles : around [-21,-22] and around [-23,^4]. We assume that the 
two cuts to be small so that we shall find a solution in a power series of the periods. 

Using the addition theorem for the prime form, we write the annulus formula ()3.33|) 
on the curve ()4.40p as 



1, E(m,m)E( 00, 00) 

— log — ^ — = ■ 

2 E{m, oo)E{m, 00) 2 Jm 



(4.41) 
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The meromorphic 1-form uJm-(x, on this curve is given by 



dz 



y{m) \ {aiz + a2)dz 



z — m 



viz) 



+ 



(4.42) 



The coefficients ai, 02 for the holomorphic forms are determined by the zero A-period 
conditions 



i=i -^^^ 



0. 



A2 



The integration in the ffist condition can be performed by deforming the contour to 
^00 U {—Cm)- As a resuh, we obtain ai = 1. The integration in the second condition is 
rather difficuh. To perform this integration, we introduce the parameterization [21] 

Zl = A21, Z2 = hA2i, 



Z3 



2 

Q + I 



— ^43; ^4 — r ^43? 



Q Zl + Z2 + Z3 + Z4 J ^3 + Z4 Zi + Z2 

^ '~ 2 ' 2 2 ■ 

By expanding the condition J^^ ujm-6b = in terms of A21, we determine the coefficient 
a2 perturbatively up to 0(A2i) as 



0(2 



^{m' - z^){m' - Z4) + 



Ail {^3 + Zi) -m' + J (m' - Z3) {m' - z^) 



^21 



AzsZiTTl' 

m' 2 (^m' + ^{m' - Z3){m' - z^)^ {Z3 + z^) {9zl - 10^3^4 + Qzj) 

+ 0{Al,), 



+ ^J{m' - Z3){m' - Z4)(^8zlzl{z3 + Z4) + 2z3ZiTn'{3zl + 2Z3Z4, + 82:4) j 



where m' := m + 



The periods of the curve ()4.4U|) 

'■^4 dz 



f-^ dz dz 

Si = 9 Tryy^)^ 82 = g — 



y{z) 



(4.43) 
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are evaluated in terms of the variables (Q, /, A21, A43). We can iteratively solve the inverse 
relations A2i(5'i, 5*2, 5-, A), A43(S'i, 6*2, 5-, A) as 6J 

A43(^i, ^2, A)2 = -^Si + -^iriC^^l - ^^^^^) + -4]^Si{hSi - 1352)(45i - 3^2) 
gA A* A' 

A2i(^i, ^2, 9, A)2 = A43(Si, ^2, 9, -A). (4.44) 

The integration of meromorphic 1-form 1)4.411) can be performed by expanding in terms 
of A21 and A43. By plugging the inverse relations ()4.44|) . we obtain the perturbative 
expansion 

-.(2) _ Si S2 







2^Am2 2^A(m + A)2 

5A2 - 8Am + Sm^ SA^ + 8Am + Sm^ „ 
^ VAW ' ^ Ag^A\m + A)^ ^ 

, 2A^ + 3A3.n - SAW - 16Arn3-8rn^^^^^ ^ ^^^3^^ 



5(2A4m3(m + A)3 

These terms precisely coincide with ()4.37p . ()4.38p and ()4.39|) . The absence of the contri- 
bution from uj^'^\z) can also be checked to all orders by numerical calculations. 



5 CFT Techniques 

So far we have considered the planar contributions to the Tri? A R corrections in super- 
symmetric gauge theories. On the other hand, there are also corrections coming from the 
torus diagrams [^3 IHH HOI HI]; which are known to be elegantly computed using CFT 
techniques ESI HI- In this section we will see how the full Tri? A R correction J-'x=o 
(= the planar gravitational corrections JF(p') ()3.34|) + the torus contribution J^f'^) is re- 
produced in this framework for the gauge theories with matter. We will first recall how it 
works in the case without matter [lH HH} HHl ISS] 1 and then examine how the matter fields 
fit in the story. 

5.1 CFT Techniques without Matter 

Consider the matrix model partition function 

Z = /rf$exp ( -— TrVr(<l>) ) (5.1) 
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with an arbitrary tree level potential 

oo 

Wi<^) = -Y.tk<^\ (5.2) 

CFT techniques are based on the the equivalence of the loop equation and the Virasoro 
constraints ISHSEHIESI- That is, defining the loop operator 

uj(z) := ^Tt^— (5.3) 



and the corrective field 



the loop equation 



^(z) := -W{z) + 2gsTTlog{z -<!>), (5.4) 
dip{z) = 2Suj{z) - W'{z), 



can be written in the equivalent form 



1 ^ ,5,,) 



\ 27ri Jc w — z ^ 

where < ■ ■ • > denotes the expectation value. Here the contour C encircles all w = A,, 
the eigenvalues of $, but not w = z. Since 

CO "1 (9 

n=0 n=-oo ^'^n 

inside a correlator, (f{z) may be regarded as a free chiral boson of conformal field theory 
and the equation ()5.5|) can be written as the Virasoro constraints 

dw^^) = 0, (5.7) 
^ c w — z I 

where 

1 
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is the energy-momentum tensor. Note that the prefactor ^ has been chosen so that their 
moments generate the Virasoro algebra in the standard normahzation. 

To leading order in ^f^, Z is computed by the saddle point approximation in which 
dip{z) is replaced with its large N expectation value —y'^^^z) T|. We need to compute 
the next-to-leading order, which may be obtained by the Gaussian approximation around 
a classical solution , and hence is described by a free conformal field theory. 

For the case of pure gauge theories without matter, CFT techniques were utihzed IS] 
to compute the genus one contribution to the Tri? A R correction. This goes as follows: 
Consider the correlation function of 2n twist operators cr(aj, aj) (z = 1, . . . , 2n) on a sphere 

gHl Uni Ea ED E21 




2n 



detA| ^JJ^Iaj — ttjl 4 ^ exp 



^m{p ■ T ■ p — p ■ T ■ p) 



(5.9) 



where r is the period matrix of the double cover with genus n (not n— 1!; the extra handle 
arises by the plumbing fixture connecting z = oo and 6b. See Section EISl) and 

Aij := <f dz^—. (5.10) 



'A. y 

The loop momenta {pi}j=i,...,n-i run over a certain momentum winding lattice. 

It was argued in , and was recently confirmed in |1U] , that the genus one free energy 
can be built from the chiral determinant of a c = 1 free boson CFT, or the chiral piece of 
()5.9|) with particular loop momenta 

n^K)) = e^^^-^detAo^ 

2n 

= e"^-"-^detA-^n(«i-«j)"^- (5-11) 

i<j 

The important point is that the twist operator a{z) does not satisfy the Virasoro constraint 
(because, for example, it has nonzero conformal weight ^) and hence must be replaced 
with an appropriately modified 'star' operator S{z) 53^ given by 

S{z) = ix._{zr^^a{z), (5.12) 
<d^{z)> = f'riai) {z - aiY-^ (5.13) 
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to this order of the approximation. If the number of cuts n is equal to the degree of 
W'{^), we find^ 

2n 

^log/i|(ai) = \ogY[{ai-aj), (5.14) 

i=l i<j 

The replacement fl5.12j) then correctly yields the genus one free energy 



2n 



, CFT 



\k=l I pi=-y/2N, 



detA ^Yliai-ttj) ^ (5.15) 



e 



as was confirmed in [211 EHl-^ Note that the CFT result ()5.15|) also automatically in- 
cludes the planar gravitational correction from sphere diagrams |16^ with the momentum 
identification^ 

Pi = -V2Ni (z = l,...,n). (5.16) 

5.2 The Full TrR A R Correction and Fermion Correlators 

We next consider the CFT interpretation of the gravitational corrections for gauge theories 
with matter. The partition function is now 



d$exp (^-^TiWi^ij e-I^^'=iT'^'°s(*- 



(5.17) 



Treating the second factor perturbatively, it is an expectation value of e~Sj=i"'"'''°s*^*~™'^) 
with respect to the measure without matter. Since 

Trlog(<l>-m7) = -^(^(m/) + iy(m7)) + Trlog(-l), (5.18) 



^If the number of cuts n is smaller than the degree of W'{^), (|5.14|l has additional contributions 
corresponding to the moment factors in eq.(4.9) of ^40', which reduce to constants in the maximally 
broken case. 

''detA reduces to that of genus n — 1 up to an irrelevant multiplicative constant in the |Ao| — s- oo limit. 
^There are reasons why the assignment (I5.16f) is natural. First, T{z) :— (Ttj^'^ has a total residue 

A^^ from matter poles [2| associated with the total ^ momentum of vertex operators (See Section IST^ ). 
whereas the difference of the residues at z = oo and 6b is 2N. Another related reason is that if Nf — 2N, 
the total momentum vanishes at infinity with this assignment, naturally reflecting the conformal bound 
of the gauge theory. 
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it should be obtained as a correlator of vertex operators charge If we take the 

convention that 

< ip{z)(p{w) >^^^ ~ -\og{z-w), (5.19) 

the prefactor of ()5.8|) must be — | so that g"^ = This motivates us to consider the 
CFT correlator 

/ 2" Nf 1 Nf \^^'^ 

I n ^i^k, -au) n e^^^^™'^ ■ e-^^^(°°) ) . (5.20) 
\fc=i 1=1 I 

This is equivalent to a correlation function of conjugate pairs of the vertex operators on 
a hyper-elliptic Riemann surface without twist operators, except that their zero-mode 
contributions differ by a factor of two In fact, to correctly reproduce the annulus 

amplitude we also need to take the normal ordering of the vertex operators on the same 
sheets : 



Nf N Nf 



CFT 



\ 1=1 1=1 I 

= (detAo)-^ E (5.21) 

where Aq^'^ is the holomorphic block with loop momenta {pi] {i = 1, . . . , n). Again, let 
us focus on a particular block of ()5.21|) with pi = —\/2Ni [i = 1, . . . ,n). Taking the 
logarithm of this block, we have 



_ / n n Nf \ 

log aJ-^'^'^ = 2{7ftJ2 njN^N, +2mJ2N,J2 '^^] 

V i,j=i i=i i=i'''^i ) 

I ^ V E{mi,6o)E{oo,mj) 

"^2,5^, ^E(m,,m^)E(oo,c5o)- ^ ' ^ 

Tij is the period matrix of the Riemann surface on which the CFT is defined (See Sec- 
tion ESI)- The first and second terms are twice of the planar gravitational corrections 
()3.35|)()3.36p . As we mentioned, the zero- mode part of the twist correlator is half of that 
on its double cover (eqs.(4.3) and (5.21) in |S21)- The last term is equal to the annulus 
amplitude J^^^ fl3.33p : by definition of normal ordering, it does not contain the correla- 
tions between pairs of points on the same sheet. The 'star'ization of the matter vertices 
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is not necessary, since the contour C in (|5.7|) does not encircle z = mi. Thus we have 
shown that the chiral block ()5.22j) together with the corrected chiral determinant precisely 
reproduces all the contributions to the full TiR A R correction for the gauge theories with 
matter. 

Remarkably, the correlator ()5.20j) closely resembles that in the bosonization of fermions 
on Riemann surfaces [SI]. This is in accord with the fact that the correlators of det($— m/) 
can be described as free fermion insertions in the soliton theory and in some special 
case they are realized as CFT correlators [Sn]- However, a significant difference here is 
that they are all for the full matrix model correlators, whereas our result concerns only 
the contributions of 0{1/N'^). Another important point is that the normal ordering is 
required in 1)5.211) . 



5.3 Non-compact Calabi-Yau as a Pinched Riemann Surface 

The A/" = 1 supersymmetric gauge theories are realized in type IIB superstring theory 
on non-compact Calabi-Yau three-folds. If one computes the effective superpotential by 
using the geometric transition, one needs to introduce a cutoff Aq to regularize some 
periods of a reduced non-compact Riemann surface. Although the quantities computed 
in matrix models are finite (since the eigenvalues are distributed on finite intervals), a 
cutoff is needed again if the free energy (or the disk amplitude) is written as a contour 
integral around infinity. Thus a matrix model curve with n cuts is a genus n — 1 punctured 
Riemann surface. 

On the other hand, the CFT correlator formula obtained in jH2] and used above as- 
sumes that the Riemann surface on which the CFT is defined is compact with no punctures. 
How can we understand the cutoff Aq here? In fact, the Riemann surface for the CFT 
can be thought of as a pinched Riemann surface of genus n, where |Ao/A|^^ is identified 
as the the pinching parameter t of the plumbing fixture. Indeed, suppose that we identify 
two annular regions near oo, 6b 



p G on the 1st sheet 



|A| < \z 



< 



pern' 



on the 2nd sheet 



|A| < \z{p)\ < 




(5.23) 
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by gluing them with a cyhnder 



S = {{X,Y) eC^ \ XY = t} (5.24) 

through the identification 

{X,Y) = (^,^1 hrWu 

(5.25) 

with \/t = |^|. Then Bi's {i = 1, . . . ,n) become closed contours, and we can verify by 
a change of the canonical homology basis that the period matrix {— J^^dz ujj}ij=i^...,n 
coincides with Fay's period matrix formula for the pinched Riemann surface [n21 in the 
|Ao| oo limit. 



|A| 


tz 


z 




tz 




Wy 


z 



6 On-shell Condition 

6.1 The M = 2 Limit and the Torus Contribution 

The M = 2 supersymmetry can be restored for U{N) gauge theories if the superpotential 
is in the form 

W^(0) = ^7A^+i (^0^^^ + E (6.1) 

so that the gauge breaking pattern is U{N) ?7(1)^; an A/" = 2 theory can be realized 
by taking the Qn+i limit [23]. The Seiberg-Witten curve is recovered from the matrix 
model curve in this limit if the moduli 5*^ satisfy the on-shell condition dWes{Si)/dSi = 
[T71ITH]. (See also Turning off the superpotential term in this way, the Af = I Weyl 

superfield Gap-^iO) and the graviphoton superfield FapiO) are combined into a single Af = 2 
Weyl superfield Wai3{0,ilj) = Fa/3{9) + ip^Gap-^. Consequently, the N = \ gravitational 
F-terms yield to an A/" = 2 F-term 

j d^xd^eT^=^{S,)W^pW^. (6.2) 

One can thus recover the M = 2 results by simply replacing the matrix model curve with 
the Seiberg-Witten curve. 
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The W'^ factor contains the Tri? A R component; it can be calculated by topologically 
twisted Af = 2 Yang-Mills theories [T3j, since the effect of the topological twist becomes 
invisible on a hyper-Kahler four-manifold M4 and the coefficient of TrR A R term of 
physical M = 2 theories agrees with that of topological theories. The latter yields up to 
a constant factor 

T^=, = b{S,) - ^c(^,), (6.3) 

1 / z^~^ 
b{Si) = --logdetAij, Aij = f dz , 

2 JAi y 

1 2N 

c{S.)= log A, A = l[{z,-z,y, 
J-^ i<j 

where b{Si) and c(S'j) are defined on the Seiberg-Witten curve 

Nf 2N 

= PM{zf - A^^-^^ - ^i) = - ^^)- (6-4) 
1=1 1=1 

On the Seiberg-Witten curve, the torus free energy which is obtained by conformal 
field theory techniques exactly coincides with ()6.3p . Recently, it was confirmed that the 
topological partition function coincides with the multi-instanton counting formula ^IJ for 
g = 1 ^31 • Thus we find consistency of our formula with the M = 2 SYM results. 

6.2 The Origin of the Planar Contributions? 

We have, however, also the planar gravitational corrections, which have no counterpart 
in the instanton calculations oi M = 2 theories. Where do these terms come from? A 
suggestive observation is that all the planar terms can be combined into a simple integral 



Y i,j=i i=i /=i-^"^ / ^ i,j=i 



E{mf, 6o)E{oo, mj) 
E{mi, mj)E{oo, 60) 



E^J. -E / ]dzT{z), (6.5) 



i=i " 1=1 



where 



I Nf 

T{z) dz = ^NiUJi + — ^Urni-^ (6.6) 
i=l 1=1 
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is nothing but the gauge theory expectation value of (Tij^j and solves the Konishi 
anomaly equations ^Hl! Since the B-period integrals of T{z) on shell may naturally be 
regarded as NS-NS fluxes on a deformed Calabi-Yau three-fold 



W'iz)^ + fn-iiz) + y'' + v^ + w'' = 0, (z, y, V, w) G C^ (6.7) 

it suggests that JF^p^) would correspond to some interaction term on the D5-branes couple 
to 2- form field B2 = Br + tBns before the geometric transition. This motivates us to 
consider the induced Chern-Simons term on Dp-branes 



/ = / C Ach(F) A \/a{R), (6. 



where C is a form field on the Dp-branes. ch(F) and A{R) are defined for the Chan-Paton 
gauge and tangent bundles on Mp+i, respectively. In our case D5-branes are wrapped 
around 2-cycles in a Calabi-Yau three- fold and filling the four-dimensional space-time. 
From this Chern-Simons term we obtain the TrR A R term 

/ C A (trl)Tri? A R, (6.9) 

JMe 

where trl comes from ch(F) and C is restricted to a 2-form in the Calabi-Yau direction. 
In this brane setup we will have the TrR A R term 

n Nf -\ 

T.^^ L +J2 ]Bm TtRAR. (6.10) 

Note that the weight factors A^j naturally arise from trl. Thus, we conjecture that the on- 
shell planar gravitational corrections may have their origin in the Chern-Simons coupling 
induced on the D5-branes. 

Although these arguments are not conclusive, the relation ()6.5|) is certainly suggestive 
and would be worth being studied. 



7 Summary and Discussion 

In this paper we have studied the gravitational corrections to the effective superpotential 
in four-dimensional A/" = 1 U{N) gauge theories with flavors in the fundamental repre- 
sentation, using the matrix model approach of Dijkgraaf and Vafa. We derived a compact 



33 



formula for the annulus contribution to the corrections in terms of the prime form on the 
matrix model curve. We also showed that the full TtRAR correction containing the torus 
as well as all the planar contributions can be reproduced as a special momentum sector of 
a single c = 1 CFT correlator. The M = 2 limit of the torus contribution agrees with the 
answer of the multi-instanton calculations and also with the geometric engineering argu- 
ment from the topological A-model. The planar contributions, on the other hand, have 
no counterpart in the instanton calculations of A/" = 2 gauge theories, and we speculated 
that the latter might correspond to the Chern-Simons term induced on the D5-branes. 

The CFT correlator we found is very close to the fermion correlator on the Riemann 
surface. In P^] the correspondence between the non-compact B-branes and fermions has 
been discussed. Since the matter contributions come from the Nj D5-branes wrapped 
around non-compact 2-cycles, our result is consistent with the picture (22 122] • In the 
topological B-model, the genus one partition function is given by a generalized holomor- 
phic Ray-Singer torsion on the Calabi-Yau geometry j2H] 



where Ap^^ is the Laplacian acting on (p, g)-forms. By using Quillen's anomaly, we also 
find the one-loop open topological string partition function. Our formula may lead to the 
explicit expression for it on a non-compact Calabi-Yau three-fold. 

There are various interesting directions to extend our analysis carried out in this paper. 
A simple extension of our analysis will be to investigate the gravitational corrections 
to the SO/Sp gauge theories. If the Af = 2 theory is deformed by an adjoint chiral 
superfield, the Chern-Simons coupling term is also induced by the orientifold plane 
It would be interesting to see whether such terms arise from the Klein bottle amplitudes 
in orthogonal matrix models. On the other hand, real symmetric/symplectic matrix 
models describe non-perturbative aspects of the SO/Sp gauge theories with flavors in 
the symmetric/antisymmetric tensor representations, respectively. The CFT description 
of these matrix models can be given by a c = —2 CFT [60j. We will report on this 
elsewhere. 

Another interesting issue to consider concerns higher dimensional gauge theories. It 
was proposed in that some five dimensional gauge theories can be described by 

unitary matrix models. The associated matrix model curve is represented by a pair of 
cylinders, where z = +oo and — oo cannot be distinguished. As a result there exist two 



3 





p,q=l 
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chiral boson fields in the CFT description. It is interesting to investigate whether the five 
dimensional multi-instanton calculations can be reproduced in this CFT analysis. We 
also expect that the meaning of the five dimensional Chern-Simons term ^1 |H2] will be 
clarified. 
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Appendix A 

In this Appendix we give an alternative proof of ()3.17|) . See Section for the definitions 

of Cj(2)(^), y^i^z), etc. 

A.l Singularities of u^'^\z)dz 

We will first determine uj^'^\z)dz from its singularities. In fact, to prove the assertion, the 
precise form is not needed. We will nevertheless derive it for future convenience for it is 
relevant to the computation of J^'^\ 

First of all, uj^'^\z)dz (as well as other higher uj^^\z)dz''?, for all k > 1) has zero A- 

periods, for the same reason as before Sj's are independent of Qs- In fact, uj^'^\z)dz 

has only singularities at the 2n branch points z = aj(0) {i = l,...,2n), which are of 
second order. They arise since the matter insertions change the locations of the cuts. We 
will show that uj^'^\z)dz is an abelian differential of the second kind with zero A-periods 
and given by 
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for some pj (j = 0, . . . , n — 2), where a, := ai(0), a'- := a[{0) and 



.(1) 1 r . w'iz) 



Mr = TT-f dz , y (z = l,...,2n). (A.2) 

m/^'*'s are the moments introduced in jHT]. p^-'s (j = 0, . . . , n — 2) are determined so that 
uj^'^\z) has zero A-periods. 

Instead of studying the singularities of uj{z, Qs) directly, it is easier to investigate those 
oiy{z,gs) ()3.29|) . y^'^\z) and —2Su;^'^\z) have the same singularities since the difference 
between the two is only linear in Qs- 

It is easy to see that y^'^\z)dz is regular at z = rrii, fhi {1 = 1,..., Nj) since J2i ^-mi 
has no singularity of 0{gl). The integral representation ()3.29p shows that y^'^\z)dz has 
possible singularities a.t z = oo, 6b and the branch points z = aj{0) {j = l,...,2n). 
However, using the conditions (j3.22|) as well as those obtained by acting ^|g^=o and 
■§^\g^=o, it can be shown that the singularities at z = oo, oo precisely cancel. Therefore, 
we have only to examine the singular behavior of y^'^\z) at the branch points. 

Since 



(2)/ A 1 d'^y{z,gs 



2 dgl 



(A.3) 



we see that the singularities at the branch points arise when ^ 



acts on the first 



factor \/n?=i(-2^ ~ o.ii9s)) of y{z,gs), and therefore y^'^\z) is proportional to 



2711 \Jc JcJ rr-i/.v/nS" 



z - w)^Ui=iiw - ai{gs)) 

1 / , W'(w) 
(p aw- 



gs=o 



-Mf^ (j = l,...,2n). (A.4) 

The contour C surrounds all the cuts, and Cz surrounds w = z. Matching the Laurent 
coefficients at the branch points, we obtain 

1 a'^M^^^ ^" 
y'^'^Hz) = = — - TT (oj — ak) + holomorphic differentials. 

(A.5) 
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This implies the equation ()A.1|) . 

Unhke to^^\z), lu^'^^z) depends on the potential not only through the positions 

of the cut but also through the moments Mj;^\ reflecting the non- universality of Tq^^ 

A. 2 An Alternative Proof of ( IXTTb 

We will now prove ()3.17|) . Using the saddle point equation 



dX J-oo \ — X 

the large N free energy (|3.25p is written as 

n / /-| | \ \ \ 

J'oigs) = Y.S,i^J_^dXp{X,gs) i^--W{\,g,) + ^log ' - -W{Xo„gs)j ■ (A.7) 

Note that ()A.6|1 is true only if A belongs to the nonzero support of p{\,gs), that is, only 
if A is such that p{\,gs) 7^ 0. Aqj is an arbitrary point on the j-th cut, whose location is 
independent of gs- 

We can extract the annulus contribution 

= (A.8) 

S - 



:= -- r^;Ap«(A)5:iog(A-m,), 

A J— 00 

^f(II) := -5:^y^rfAp(2)(A)(H^(A)-2^1og|A-Ao,|). (A.9) 

j=i °° 

(log A drops out since dXp'''^\X) = 0.) Since J^q'\i) is already the right hand side of 
(|3.17|) . we must show that J-'q^\iI) = 0. For this purpose let us first rewrite ()A.9|) in a 



contour integral on the first sheet using uj^'^>{z). We take the branch cut of log(2; — A 
so that it runs from 2; = A + iO (2; = A on the upper side of the cut)to z = —00. Then 
depending on z < j, = j or > j we have 

rfAp(2)(A)log|A- Aojl 



1 P 

— dz J^\z) {\og{z - Xoj) + Hi) {i<j), 
zm J Ai 

7^ i f dz J^\z) \og{z - Ao,) + m /'"'~'° dz J^\z)] {t = j), (A.IO) 

— [ dzuj^^\z)log{z- Xoj) {i>j), 
y Zm JAi 
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where the contour Aj(Aoj) starts from z = Xqj + iO, surrounds the jth cut once anti- 
clockwise and goes back to 2; = Xqj + iO again but on the other side of the branch cut of 
log(z — Xoj). From ()A.10|) we find 



/oo 1 /■ 

rfAp(2)(A)log|A-Ao,| = -^,lim I dz y^'\z) (j = l,...,n). (A.ll) 
-00 4o \Ao\^ooJBj 

Therefore, to prove Fq^'*(II) = it suffices to show that 

1 n „ n „ 

J2f dz/'\z)Wiz)+Y.sJ dzy^'\z) = 0. (A.12) 

This can be proved by applying Riemann's bilinear identity, as we will show below. 
Cut out along 

A„_i5„_iA;iifi-ii ■ ■ ■ A,B,A^'B^\ (A.13) 

Sn-i is represented by a 4(n — l)-sided polygon with identification in a standard manner. 
Since y^'^\z)dz is an abelian differential of the second kind, its integral 

Y{z) := r dxy^^\x) (A.14) 

JzQ 

defines a single-valued meromorphic function Y{z) inside the polygon for an arbitrary 
reference point zq. Let us evaluate the integral of 

T] ■.= Y{z)y^^\z)dz (A.15) 

along the 4(n — 1) sides of the polygon. Applying Riemann's bilinear identity, we obtain 

2m V ReSpT] = V f / dzy^'^\z) l dxy^^\x) - i dzy^'^\z) i dxy'^'^\x)] . 
p^s^s_^ j=i V^j Jb, Jb, J a, J 

(A.16) 



The right hand side is 

n-l 

IB, 



n—1 „ 

right hand side of (|Xl6l) = Am ^ Sj f dz y^^\z). (A. 17) 
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Let us compute the left hand side. Since 



9 

z 

,^00 _^^_^/^^)^(^(^-2)^ (A.18) 

we find 

27riReSoo?? = / + 47r2^r (00), 

27riRes6b?? = -/ - 47rz5r (60) 

= +i dzy^^\z)W{z) -AmSY{6o), (A.19) 

where is the same as A^o but on the second sheet. We have shown in section lA.ll 
that y^'^\z)dz has second order poles at every branch point z = ai (i = 1, . . . ,2n), and 
therefore Y(z) has first order poles there, while we also see from ()H.29j) that y^^^dz has 
first order zeroes at those branch points. Therefore r] is regular there and hence has no 
nonzero residues. Thus, in all, we have 

left hand side of ^K^^ = 2<f dz y^^\z)W{z) + AiriS {Y{oo) - Y{6o)) . (A.20) 
By definition 

POO 

Y{oo)-Y{6o) = / dzy^^\z), (A.21) 



where the contour cannot cross any side of the polygonal region since Y{z) is single- valued 
only inside. Therefore the contour can pass only through the nth cut: 

r(oo) -r(6b) = lim / dzy^^\z). (A.22) 

Plugging ()A.22|) into ()A.20|) and equating it with (jA.17|) . we obtain the desired equation 
(HHl). 

Appendix B The Special Geometry Relation 

It is now well-known that the sphere amplitude J-'jf^ satisfies the special geometry relation 



hm (- i dzy^'\z)-W{k,)+2S\og—A . (B.l) 



dSi lAoHoo V2 Jb, 
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This fact was pointed out and proven using the energy cost argument in jT] . A proof using 
a Legendre transformation was given in [^. In this Appendix we will give a sketch of 
an alternative proof neither using Legendre transformations nor resorting to any physical 
argument. 

Before we consider (IB. II). we first recall the relation 



-dz 



-AmCji{z) (i = 1, . . . , n), 



(B.2) 



where, as we defined in the text, y^^\z) is related to the large N resolvent oj'^^^z) by 

= W\z)~2Suj^''\z). (B.3) 

Note that the relation ()B.2|1 is also true when the degree of W'{z) exceeds the number of 
the cut n and y'^'^\z) has extra zero factors (See footnote |21). To see this we first use the 
integral representation of y'^^\z) (obtained by setting (7^ = in ()3.29|) ) to verify that the 
possible singularities of ^^gg*-^-* are only located at z = 00 and 6b. Next we write 



yW(^) =^w'{zy + f{z,{s,}) 

for some /(z, {Sj}) to find that 



W'(z) 



Oil) 



[B.4) 



(B.5) 



smce 



= W'{z)-^W'izr + fiz,{S,}) 

w'{z) + ^w'{zy + f{z,{Sj}) 

must behave like — + 0{^) as 2 ^ 00. Therefore we find that 



:b.6) 



dy^ 



dz 



2JW'izr + fiz,{S,}) 



-dz 



:b.7) 



has at most a first order pole at 2; = 00, and hence can be written as a linear combination 
of cDj's. Then flB.2|) follows from the fact that the Aj-period of y^^\z)dz is —AniSi (i = 
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We will now prove the special geometry relation ()B.1|) . ^o°^ is given by 

S /•« 



(0) 







rfAp(°)(A)lo^ 



lA-A, 



Oil 



A 



(B.8) 
(B.9) 



We express it in terms of contour integrals of the complex plane. Note that the integral 
of dzy^^\z) \og{z — Xoj) along A^o does not vanish as |Ao| — > oo, and hence must be taken 
into account: 



(0) 



lim 

lAoHoo 



Siri J A 



dzy^^\z)W{z) 



n Q 



A 4 Jb. 

Differentiating it by Si, we find 



dz y^'\z)-^W{Ao) + SHog^ 



(B.IO) 



(0) 



lim 

|AoHoo 



I r n 
-f u,{z)W{z)-7czY.S, I 

+ \ dz y^'\z) - ^W^(Ao) + Slog ^ 



a),(z) + S'log^^ 

Bi A 



(B.ll) 



where the S'logAo terms are separated so that each line gives a finite result. Again, we 
can simplify it by using Riemann's bilinear relation. In this case we cut out the Riemann 
surface along 

An-iBn-iA-\B-^, ■ ■ ■ AiBiA^'B^'A^BnA^iBn)-' (B.12) 

and evaluate 

dz y''^\z)Qi{z) , Qi{z) := f uJi{x) (i = l,...,n), uJn{z) := uJn{z) (B.13) 
along the boundary of the resulting polygon. The result we obtain is 



lim 

|Ao|^oo 

for ? = !,.. 



/ i2;,{z)W{z) + 2m Sj ! uj,{z) + ]: 1 dz y^''\z) - W^(Ao) 
Ja^ jr[ JBj 2 JBi 



(B.14) 



, n. 



Using ()B.14|) in ()B.11|) . we obtain the special geometry relation ()B.l^ 
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Appendix C The Large N Free Energy for a Quadratic 
Potential 



^0 



where 



S 
2 



S 



gat 
s 



a{gs) - a{gs) ^ 



gstfi{gsY A (l^{gs)a{gs 



+ 



2S V'"°V 2 ; ' 2a((7,)2 

gj I a{gs)-\\og{-i^) - 1) + a{gs) \oga{g,) 
' S 



a{9s) - a{9s) ^ 



- log(a(fi'.) - a{gs) ) 



+tfi{gs)^ log 



K9s) , c{gsf 



K9s 



+2^ 



+ 



ai9s) - a{gs) ^ 
+ 9s log(-m) 



1 m 



(C.l) 



a{9s 



m - c{gs) ± ^/(m - c(5(,))2 - fi{g,y 
K9s) 



(C.2) 
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